Stiffness of Contacts Between Rough Surfaces 
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The effect of self-affine roughness on solid contact is examined with molecular dynamics and 
continuum calculations. The contact area and normal and lateral stiffnesses rise linearly with the 
applied load, and the load rises exponentially with decreasing separation between surfaces. Results 
for a wide range of roughnesses, system sizes and Poisson ratios can be collapsed using Persson's 
contact theory for continuous elastic media. The atomic scale response at the interface between 
solids has little affect on the area or normal stiffness, but can greatly reduce the lateral stiffness. 
The scaling of this effect with system size and roughness is discussed. 

PACS numbers: 46.55.+d, 62.20.Qp, 81.40.Pq 



The presence of roughness on a wide range of length 
scales has profound effects on contact and friction be- 
tween experimental surfaces. Under a broad range of 
conditions [3-@, the area of intimate contact between 
rough surfaces A c is orders of magnitude smaller than 
the apparent surface area Aq. As discussed below, this 
provides the most common explanation for Amontons' 
laws that friction is proportional to load and indepen- 
dent of Aq. Because A c is small, the interfacial region 
is very compliant. In a range of applications the inter- 
facial compliance can significantly reduce the stiffness of 
macroscopic joints formed by holding two components 
together under pressure [l], 0j . 

In this paper, we examine the effect of surface rough- 
ness on the normal and lateral stiffness of contacts be- 
tween elastic solids using molecular dynamics (MD) and 
continuum calculations. The results provide a numeri- 
cal test of recent continuum theories [H, Q and their ap- 
plicability to real solids. The contact area and normal 
stiffness approach continuum predictions rapidly as sys- 
tem size increases. Continuum theory also captures the 
internal deformations in the solid under tangential load, 
but the total lateral stiffness may be greatly reduced by 
atomic scale displacements between contacting atoms on 
the opposing surfaces. This makes it a sensitive probe 
of the forces underlying friction and may help to explain 
unexpectedly small experimental results [Toj | . 

The topography of many surfaces can be described as 
a self-affine fractal 0, El- Over a wide range of lengths, 
the root mean squared (rms) change in height dh over a 
lateral distance t scales as a power law: dh ~ £ H , where 
the roughness or Hurst exponent H is typically between 
0.5 and 0.9. Greenwood and Williamson (GW) consid- 
ered the peaks of rough landscapes as independent asper- 
ities and found that A c rose linearly with normal load Fjq 
for nonadhesive surfaces Q. This explains Amontons's 
laws if there is a constant shear stress at the interface. 
A linear scaling of area with load is also obtained from 
Persson's scaling theory, which includes elastic coupling 
between contacts approximately 12, 13| . 



Dimensional analysis implies that the linear relation 
between load and area must have the form 



A C E' = kF n /^/WW 



(1) 



where a modulus like the contact modulus E' is the only 
dimensional quantity characterizing the elastic response, 
and the rms slope the only dimensionless quantity charac- 
terizing the roughness. Numerical solutions of the contin- 
uum equations [4, 6] show that n is near 2. Results for dif- 
ferent H and Poisson ratio v lie between the analytic pre- 
dictions of GW, \/2n ~ 2.5, and Persson, ~ 1.6. 
One advantage of Persson's model is that, as in numeri- 
cal results, A c /Fn is constant over a much larger range 
of loads than GW 14]. Another is that it captures [§] 
the power law scaling of correlations in contact and stress 
that was found in numerical studies 15, 16j |. 

The normal stiffness is related to the change in average 
surface separation u with load. Experiments 17|, [18[ and 
calculations 0, [H, [lj| show an exponential rise in load 
with decreasing u, Fpf = cA E' cxp[— u/jh Yms ], where 
h lms is the root mean squared (rms) variation in surface 
height and 7 a constant of order 1. Differentiating leads 
to an expression for the normal interfacial stiffness: 



k T N = —dF^/dii — Fn /7/ir 



(2) 



For self-affine surfaces this interfacial stiffness decreases 
as h~,£ a ~ L~ H with increasing system size L. Our sim- 
ulations test this scaling and show that 7 is nearly con- 
stant. They also examine the connection between this 
normal stiffness and the transverse stiffness k T at forces 
lower than the static friction [201 ] . 

We consider nonadhesive contact of a rigid rough solid 
and a flat elastic substrate. This can be mapped to con- 
tact of two rough, elastic solids in continuum theories 
0, E2] • The mapping is only approximate for atomic sys- 
tems 2l|, |22[, but working with one rigid solid reduces 



the parameter space in this initial study. Rigid surface 
atoms are placed on the sites of a simple cubic lattice 
with spacing d' , and only atoms in the outer layer are 
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close enough to interact with the substrate. They are 
displaced vertically to coincide with a self-affine fractal 
surface of the desired H . Surfaces with roughness on 
wavelengths from Z m ; n to l max were generated as in Ref. 
[I5J]. The rms slope yJ\Vh\ 2 — 0.1 for the results shown. 
Consistent results were obtained for slopes from 0.05 to 
0.15. Slopes of 0.2 or greater led to plastic deformation 
in the substrate. Continuum calculations also show plas- 
ticity for large slopes 0. 

The elastic substrate is an fee crystal with nearest- 
neighbor spacing d and a (100) surface. Surface atoms 
form a square lattice that is rotated by 45° relative to 
that of the rigid surface. This rotation and the choice 
of did! prevent commensurate locking of the surfaces 
2^, 24 1 . Substrate atoms separated by r interact with 
a Lennard- Jones (LJ) potential: Ulj = ^e[{a/r) 12 — 
(<r/r) 6 ], where e and a are the bonding energy and diam- 
eter. To speed calculations, the potential and force are 
interpolated smoothly to zero at r c = 1.8er (25)]. Since we 
compare to continuum theories with no interfacial adhe- 
sion, purely repulsive interactions are used between sub- 
strate and rigid atoms. An LJ potential with length a' is 
truncated at the energy minimum, 2 1 / 6 er'. Unless noted 
a' = 0.5cr and d! = 0.3d to minimize interfacial friction 
and speed the approach to continuum theory Q . 

Substrate atoms are arranged to form a cube of side L. 
Periodic boundary conditions are applied in the plane of 
the surface and bottom atoms are held fixed. The elastic 
constants C33 = 70.2e/cr 3 and C44 = 41.8e/cr 3 were mea- 
sured by displacing the top surface. Cubic crystals do 
not exhibit the isotropic elasticity assumed in continuum 
theory. To test for any effect from anisotropy, we also per- 
formed Greens function MD (GFMD) [26[ for isotropic 
continua with v — and 0.35. Since thermal fluctuations 
are ignored in continuum theory, we consider the limit of 
zero temperature T. The energy is minimized for given 
external forces or displacements. The fractional contact 
area A c /Aq is obtained from the fraction of surface atoms 
that interact with the rigid surface [2lj |. 

In all cases studied, A c rises linearly with Fn. More- 
over, the value of k approaches previous continuum re- 
sults as system size increases. The decrease in k with 
increasing L is similar to that found by Campana and 
Miiser in their GFMD calculations [6j. We found k con- 
verged more rapidly when their atomistic Greens function 
was replaced by an ideal elastic Greens function and we 
use this in the GFMD results below. As in previous 2D 
atomistic calculations 21], the stress and contact corre- 



lation functions from our LJ calculations exhibited the 
power law scaling q~( 1+H ) found in GFMD calculations 
and Persson's theory [t| [l6[ . 

Figure Q] shows the variation of F/v with interfacial 
separation u for several L and H = 0.5 and 0.8. In all 
cases, Fjv rises exponentially over a range of loads that 
corresponds to fractional contact areas between 1 and 
10%. Statistics are too poor at lower areas and non- 
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FIG. 1: Logarithm of load as a function of (no —u)/h m s, and 
linear fits corresponding to 7 = 0.48. The separation at first 
contact, uo, is shifted slightly to prevent overlap. Atomistic 
results are for H = 0.5 with L = 378.4<r (circles), 189. 2a 
(squares), and 94.6<r (triangles) and for H = 0.8 with L = 
189. 2a (crosses) and 94.6a (pluses). Filled triangles are for a 
GFMD simulation with L = 128d and v = 0. 



linear corrections to Eq. Q] are seen at larger areas [4| ■ 
The linear fits to all results have the same slope, corre- 
sponding to 7 = 0.48 and best fit values for all H and L 
studied differ by less than 10% from this value. GFMD 
results were at the higher end of this range and showed 
no change as v increased from to 0.35. Earlier contin- 
uum calculations [Ij], elastic atomic calculations [l9| and 
experiments [l8| were consistent with 7 w 0.4. This rep- 
resents a compelling success of Persson's approach, and 
raises the question of whether 7 may have a unique value 
in the thermodynamic, isotropic limit. 

The normal stiffness from Eq. [2] includes a component 
from the increase in contact area with load as well as 
the change in force at fixed area. There is also a compli- 
ance associated with changes in the separation between 
contacting rigid and substrate atoms that is generally ne- 
glected in continuum theory. To isolate the stiffness as- 
sociated with deformation within the substrate at fixed 
contact area k Is , we applied constraints directly to the 
atoms that contacted at a given load. The normal and 
lateral stiffness were then obtained by displacing these 
contacting atoms in the normal or lateral direction and 
measuring the change in force. The contribution from 
the bulk response was subtracted so that the stiffness re- 
flects the change in surface separation u or lateral surface 
translation ut- This approach is straightforward to im- 
plement in experiments and was found to be consistent 
with direct averaging of atomic separations. 

Figure (5Ja) shows the scaled normal interfacial stiff- 
ness k 1 ^ h rms / 'AqE' as a function of the dimensionless 
load Fn/A E' used to find the contacting atoms. Once 
again, results for all systems show the same behavior, 
and the stiffness rises linearly with load as predicted by 
Eq. [2] The points lie slightly above the dashed line cor- 
responding to 7 = 0.48 due to small deviations from the 
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FIG. 2: The scaled (a) normal stiffness and (b) tangential 
stiffness as a function of W/AoE*. Results are for H — 0.5 
(open symbols) and H — 0.8 (filled symbols) with L = 189. 2a 
(circles), L — 94.6a (squares) or L = 47.3a (triangles). 
Dashed lines have slope I/7 with 7 = 0.48. 



analytic form of Eq. [5] One might expect k 1 ^ to be 
substantially less than the total stiffness because it does 
not include the stiffness from increases in contact area. 
However, the two stiffnesses are nearly the same because 
newly contacting regions carry the smallest forces. 

The incremental response of an ideal elastic solid does 
not depend on any pre-existing deformation. This implies 
that we should obtain the same stiffness by displacing the 
same set of atoms on the initial undeformed surface. Di- 
rect evaluation of the stiffness in this way gave slightly 
lower stiffnesses than Fig. [2 with the difference increas- 
ing from the numerical uncertainty to about 15% with 
increasing F^. This provides an estimate of the contri- 
bution that anharmonic effects may make to the stiffness 
of real materials at the rms slope used here. 

The above results imply that the stiffness of elastic 
solids at fixed contact area is uniquely determined by 
the distribution of contacting points. This conclusion 
may seem at odds with Eq. [3J since the contact area 
has no independent connection to load or surface rough- 
ness. The resolution is that variations in load and rough- 
ness cancel. If the response is linear, one can scale h rms 
and Fn by the same factor and the contact area will be 
unchanged. Indeed one can combine Eqs. [T] and [5] to 
eliminate F^: 



k 1 



A E' v^jvhp 



K7 Aq' 



(3) 



For a self-afhne surface, the ratio h Tms / -\/|V/i| 2 cx 
GmaxAmirO^min depends only on the small and large 
scale cutoffs in roughness. 

Figure[3][a) shows the scaled stiffness k^ vs. area. The 
results were obtained by displacing atoms from their po- 
sitions on the initial flat surface to eliminate anharmonic 
effects. Results for all systems collapse onto a common 



straight line, providing clear evidence for the direct con- 
nection between stiffness and contacting area. The slope 
is near unity as expected from the separate values of k 
and 7. 

All of our atomic simulations show kj? /kjf > 1. This 
is surprising given that Mindlin [20j j and recent work 
[III predict k^/k* N = 2(1 - u)/(2 - u) < 1. How- 
ever, this result is for isotropic systems. One measure 
of the anisotropy of the LJ crystal is that the ratio 
C44/E' ss 0.57, while it is (1 — v)/2 < 1/2 for an isotropic 
solid. This is also consistent with shear stresses hav- 
ing a higher relative stiffness than expected. In gen- 
eral, the total elastic energy stored in the interface is 
T lq f{—q)G(q)f(q)/2 where G is the Greens function ma- 
trix relating displacements to forces / (26[. The ratio 
of stiffnesses can be obtained by averaging the diago- 
nal components of qG{q) corresponding to normal and 
tangential displacements over q and assuming the same 
power spectrum describes the respective forces. This ra- 
tio agrees with Mindlin's result for isotropic systems, but 
will vary with crystal anisotropy. 

As noted above, kjf and k^ only consider the stiffness 
associated with deformation inside the substrate at con- 
stant area. There is an additional stiffness k Ia associated 
with relative motion at the interface between contacting 
substrate atoms and the rigid surface. Simulations of sin- 
gle asperity contacts show that this atomic scale motion 
has very different effects on the total normal and tangen- 
tial contact stiffness [22j. The stiffness k 1 ^ resisting nor- 
mal displacements is large, because contacting atoms are 
in the steep repulsive part of the L J potential. The com- 
pliance associated with these bonds has little effect on 
the total stiffness of single asperities or the multiaspcrity 
contacts studied here. As shown above, the full normal 
stiffness including this compliance and changes in area 
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FIG. 3: The scaled (a) normal stiffness and (b) tangential 
stiffness as a function of A c /Ao. Results are for H — 0.5 (open 
symbols) and H — 0.8 (filled symbols) with L = 189. 2a (cir- 
cles), L = 94.6a (squares) or L = 47.3a (triangles). Crosses 
show GFMD results for v — 0. Dashed lines have slope 0.87. 
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(Fig. [T]) is consistent with the stiffness from compression 
of the elastic substrate alone (Fig. [2]). 

In contrast, the tangential stiffness of single asperity 
contacts can be an order of magnitude or more smaller 
than continuum predictions 22] . Surfaces rarely share a 
common period, and the force per area preventin g la teral 
motion averages to zero as the contact size grows 23|, 24 1 . 
In our simulations, k^ a depends strongly on d'/d, er'/er, 
l m m, and the orientations of the solids, and the total 
lateral stiffness is generally too small to distinguish on the 
scale of Fig. [5J Indeed, this is desirable in comparing to 
continuum theories for normal contact since they assume 
zero friction and lateral stiffness. 

Our results for a wide range of parameters can be sum- 
marized by noting that Eq. [1] arises because the distri- 
butions of pressure and contact sizes are independent of 
load. This implies k^ a = otA c c^j <j where a is of order 
one if the contact has the same stiffness as the substrate. 
Since compliances add, the total interfacial compliance is 



1 = 1/ k^ a + 1/ kip s . Both contributions diverge as Fm 
and A c go to zero, explaining why 1/fcp can be compara- 
ble to the bulk compliance of the bounding solids [l|, 0] ■ 
Since both contributions scale with A c but 

with a lower slope than in Fig. [31 

The substrate compliance dominates for suffi- 
ciently large systems and roughness since k^ a /k^f = 
ac4 4 h rms / '<jE' h\ 2 ~ (/maxAmin) H - However, for 
small a, the value of Z max required to reach this limit 
may be large, particularly on the scale of atomistic sim- 
ulations, microclcctromechanical systems, or the wave- 
length of ultrasound used to measure shear stiffness [l(| ■ 
For the case of H = 0.5 in Fig. [3j k I Ta still reduces the 
total interfacial stiffness by a factor of two for l max ~ 
1300cr ~ 400nm with a = 0.1 and for l max ~ 40/j,m for 
a = 0.01. Our measured a span this range and one can 
estimate a in experimental systems from the static fric- 
tion coefficient fj, s . If /3cr is the lateral displacement for 
the force to reach the static friction, then the above equa- 
tions yield a — fi s y/\Vh\ 2 /K0. Typical values of /? are of 
order 1/4 so a and fi are of the same order. 

In conclusion, atomic scale simulations were used to 
study contact between surfaces with roughness on a wide 
range of scales. The results for area and normal stiffness 
are consistent with Persson's continuum theory down to 
relatively small scales even though the solid is not contin- 
uous or perfectly elastic. The area and internal stiffnesses 
of systems with a range of H, L and v show the linear 
scaling predicted in Eqs. [T]-|3]with nearly constant values 
of k and 7. The internal stiffnesses were shown to depend 
only on the geometry of the contacting region. Atomic 
scale displacements between contacting atoms have little 
effect on the normal stiffness, but can change the lateral 



stiffness by orders of magnitude. This sensitivity makes 
lateral stiffness a promising probe of the atomic scale in- 
teractions that underly friction. 
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